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FORMATION  AND  ADIABATIC  COMPRESSION  OF 
REVERSED-FIELD  THETA  PINCHES  IN  IMPLODING  LINERS 

I,  Introduction 

In  the  NRL  "Captive  Liner"  concept, a liquid  metal  liner  Im- 
plodes, heating  a plasma  by  adiabatic  compression  to  thermonuclear  con- 
ditions. The  liner,  which  Is  stabilized  against  Raylelgh-Taylor  In- 
stability by  rotation,  rebounds  after  Ignition,  and  the  process  Is 
repeated  cyclically. 

One  Important  problem  under  study  Is  that  of  creating  a suitable 

plasma  In  the  Interior  of  the  device  during  each  cycle  prior  to  the 

Implosion.  Access  to  this  region  Is  possible  only  through  the  ends. 

The  plasma  must  be  sufficiently  conducting  (T^  ^ 100  eV)  that  a buffer 

magnetic  field  can  be  maintained  between  plasma  and  liner  during  the 

Implosion.  Further,  since  In  most  variants  of  the  concept  the  Implosion 

Is  quite  slow  ^ 1C  cm/sec),  the  plasma-field  configuration  must  be 

hydromagnetlcally  stable  over  times  > 10  ^ sec. 

A promising  candidate  to  meet  these  requirements  Is  a reversed  field 

theta  pinch  created  by  a rotating  relativistic  electron  beam  Injected 

through  the  end  of  the  system  (see  Fig.  1).  This  configuration  has  been 

T2l 

demonstrated  experimentally.*-  - Reversed- field  theta  pinches  (belt 
pinches)  have  displayed  exceptional  stability,  apparently  due  to  finite 
ion-gyroradlus  effects.-  - 

To  evaluate  the  prospects  for  e-beam- induced  reversed-field  con- 
figurations, several  questions  must  be  answered.  One  concerns  the 

Note;  Manuscript  submitted  May  2,  1977. 
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scaling  laws  relating  the  e-beam  parameters  to  those  of  the  configura- 
tion (temperature,  density  and  field  profiles).  Experiments  now  In 
progress  at  NKL  are  aimed  at  detemlnlng  these  laws.  A second,  the  one 
addressed  here.  Is  related  to  the  changes  In  the  equilibrium  resulting 
from  compression,  and  their  effect  on  the  system' s efficiency  as  a 
reactor.  Logically,  the  next  question  (which  has  not  yet  been  considered 
to  any  extent)  would  be  the  stability  of  such  configurations.  Presumably 
even  Initially  stable  configurations  would  become  less  so  under  com- 
pression as  the  density  and  field  gradients  steepen,  and  for  any  distri- 
bution a threshold  must  exist  beyond  which  further  compression  produces 
unacceptable  or  even  catastrophic  destabilization.  The  characteristic 
time  over  which  plasma  confinement  at  maximum  density  must  be  maintained 

Is  .V/  a /V,  , where  a Is  the  final  (minimum)  liner  radius.  This  Is  sub- 
o L’  o 

stantlally  shorter  than  the  total  Implosion  time  (~  1 msec). 

Throughout  this  paper  we  ignore  end  effects  and  azimuthal  dependence, 
so  that  field  and  fluid  quantities  depend  only  on  r.  To  give  an  idea  of 
the  sort  of  behavior  to  be  expected,  we  begin  by  considering  briefly  the 
simplest  Ideal  case,  that  of  unmixed  plasma  and  flux.  In  such  an  ideal 
final  compressed  configuration,  the  plasma,  field  and  liner  regions  are 
distinct  and  sharply  demarcated  (Fig.  2).  Pressure  balance  requires 


that 


Po  = 


2U. 


(1) 


where  p^  is  the  plasma  pressure  and  the  strength  of  the  (uniform) 
buffer  field.  The  region  occupied  by  the  reversed  field  Is  assumed 
vanishingly  small.  This  Is  only  valid  If  a relatively  small  fraction  of 


ii 


the  field  lines  of  the  total  system  are  closed.  Even  In  this  limiting 
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case,  the  magnetic  buffer  zone  must  have  a finite  extent  (b^  > a^).  Its 
thickness  must  be  at  least  several  ion  gyroradli  or  mean  free  paths 
(if  a cold  background  plasma  is  present)  in  order  that  the  fluid  picture 
be  valid,  and  it  must  be  thicker  than  the  (as  yet  undetermined)  region 
within  which  copious  vapor  or  plasma  from  the  liner  is  present.  On  the 
other  hand,  if  it  occupies  a volume  which  is  not  small  compared  with  the 
total  compression  volume,  an  Intolerably  large  fraction  of  the  system 
energy  is  diverted  from  the  resulting  plasma  payload. 

The  D-T  reaction  rate  divided  by  averagedover  a Maxwellian  dis- 
bution  has  a broad  maximum  at  about  12  keV.  Hence  the  payload  state 
should  Ideally  have  a uniform  temperature  profile  which  attains  this 
value  at  peak  compression,  as  well  as  a uniform  density.  Given  a con- 
figuration of  this  description,  it  is  immediately  possible  to  deduce  the 
starting  conditions  (prior  to  compression)  which  will  give  rise  to  it. 

Let  the  precompression  radius  of  the  liner  be  called  b^,  that  that  of 
the  plasma-field  Interface  a^.  Then  conservation  of  flux  Implies 

Bi  [b/  - ./]  (2) 

Adiabatic  compression  of  the  plasma  implies  an  initial  pressure  (assuming 
Y = and  2-dlmensional  compression)  given  by 


whence  the  magnetic  field  strength  6^  is  found  from 


conservation  of  total  particle  number  Implies  a number  density 

which  together  with  Eq.  O)  gives  a temperature 
^1  - 


These  relations  suffice  to  determine  the  eturller  state  entirely,  given 
the  final  values  a , b , n and  T (or  B or  p ),  and  one  other  quantity, 
e.g.,  the  compression  ratio  a = b^/b^.  For  example,  we  have  the  results 
shown  In  Table  I. 

Table  I 


a (cm) 
b (cm) 

a 

-5 

n(cm  ) 
T(eV) 
p (kbar) 
B (kG) 

€ 


Final 

State 

Initial  State 

1.0 

20.0 

30.0 

4o.o 

50.0 

1.5 

24.2 

55.5 

46.7 

57.9 

1.0 

16.1 

25.0 

51.2 

58.6 

2.1(-8xl0^® 

6.20x10^^ 

2.76x10^^ 

1.55x10^^ 

14 

9.92x10 

lo'^ 

184 

107 

75.1 

5i^.5 

59.5 

1.81 

0.468 

0.180 

0.08; 

10^ 

6.78 

5.45 

2.14 

1.47 

0.667 

0.765 

0.789 

0.804 

0.816 

S' 


Here  € is  the  ratio  of  plasma  energy  to  total  energy. 


6 = 


a^  ^ n k T 
2 


(7) 


I n k T + (b^-a^ 

It  is  clear  that  as  we  retreat  to  earlier  and  earlier  initial  states, 
the  buffer  region  becomes  thin  relative  to  the  plasma  radius,  although 
Its  absolute  thickness  Is  greater  then  in  the  final  state. 

The  following  phenomena  act  to  prevent  attainment  of  an  Ideal  com- 
pletely separated  final  plasma-field  configuration:  particle  and  thermal 

diffusion  out  of  the  plasma;  Inward  diffusion  of  liner  vapor  and  plasma; 
radiative  cooling  of  the  plasma;  and  Initial  deviations  from  complete 
separation  due  to  transport  processes  and  to  the  finite  annular  thickness 
of  the  Injected  e-beam.  The  last  of  these  represents  the  particular 
concern  of  the  present  report.  We  develop  a formalism  suitable  for 
determining  the  successive  stages  through  which  a mixed  plasma-field  con- 
figuration must  pass  under  adiabatic  compression,  and  apply  It  to  a simple 
but  realistic  example. 

The  plan  of  the  report  is  as  follows.  In  Section  II,  we  summarize 
the  plasma  processes  (Including  anomalous  heating)  Involved  In  producing 

the  Initial  current  distribution.  In  Section  III  we  Introduce  a 

\ 

Lagranglan  formalism  to  describe  the  compression  history  and  apply  It 
to  a configuration  derived  In  Section  II.  We  conclude  in  Section  IV 
with  a discussion  of  our  results  and  some  recommendations  for  future 
study. 
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II.  Formation  of  Field-Reversed  Equilibrium 

The  experiments  presently  underway  to  study  the  production  of  a 
reversed-fleld  configuration  using  relativistic  electron  beams  are  still 
in  the  preliminary  stages.  Therefore,  a detailed  characterization  of 
the  magnetic  field  and  plasma  pressure  profile  Is  not  yet  available. 
However,  results  to  date^^"^  cb  Indicate  that  the  beam-plasma  interaction 
Is  highly  anomalous,  since  much  higher  plasma  temperatures  and  net  mag- 
netic fields  are  observed  than  can  be  predicted  on  the  basis  of  classical 
dissipation  of  the  beam-induced  azimuthal  return  current.  For  purposes 
of  estimating  the  characteristics  of  the  Initial  plasma  that  we  might  be 
able  to  produce  using  the  Gamble  II  beam  (electron  energy  1.0  MeV,  total 
current  1.5  MA,  pulse  length  JO  nsec),^  we  have  used  the  following 
model.  A 1-5  kA/cm^  beam  is  injected  into  a partially  ionized  hydrogen 

plasma,  inducing  an  equal  and  opposite  return  current.  This  return 

2 

current  deposits  energy  in  the  plasma  at  a rate  , where  is 

the  classical  resistivity  due  to  Coulomb  and  electron-neutral  collisions, 
plus  the  resistivity  due  to  any  instabilities  which  may  be  operative. 
When  the  beam  pulse  is  over,  the  plasma  must  reach  ^ lOOeV  temperature, 
and  the  net  current  left  behind  in  the  system,  equal  to  the  net  current 
in  the  system  due  to  return  current  dissipation  at  the  end  of  the  beam 
pulse,  must  be  large  enough  to  produce  the  desired  reversed-field  plasma 
confinement  geometry. 

In  order  to  determine  the  potential  validity  of  this  model,  a 

computer  code  which  follows  the  evolution  in  time  of  plasma  density  and 

[5l 

temperature  was  used.  This  code,  described  elsewhere,  includes  the 
principal  ionization  and  radiation  processes  and  an  energy  input  term 


6 


rj  j . Here  rf  Is  determined  at  each  Instant  of  time  from  the  plasma  con- 
dition (<n  , T , T,)  for  the  various  Instabilities  assumed  to  be  present 
'e  e 1 

In  a given  run.  Since  the  lower  hybrid  drift,  Ion  acoustic,  and  electron 
electron  two-stream  Instabilities  are  considered  likely  candidate  sources 
of  the  anomalous  heating  In  the  experiments,  we  have  tried  resistivities 
derived  for  these  three  modes  In  various  combinations  and  separately 
(but  In  ail  cases  added  to  the  classical  resistivity). 

The  resistivity  values  used,  expressed  in  ohm-cm,  were 

^ /s 

= l.Oxlo"^  T^"  ^ [tn  A + l.l^xio"^  ^V"e^  ^ 

%d  = 5 (9) 


« = 1.2x10“  n 

'M  a a 


» ^ ^crit 


= 0 otherwise; 


T,_  = 3.3x10^9  j/  ^ 


where  -tn  A is  the  Coulomb  logarithm;  n^  and  n^  are  neutral  and  electron 

-3 

number  densities  In  cm  ; B Is  magnetic  field  strength  In  kG; 

*s  (2kT^/M^)^^^  is  the  ion  thermal  velocity  in  cm/ sec;  and  are 
electron  and  ion  temperatures  in  eV;  v^  Is  the  electron  drift  velocity, 


v_  = 6.0x10  J/n  cm/sec; 
D e 


Vcrit»  the  lon-acoustlc  cutoff.  Is  given  by 
\rit  “ 9.8x1o5  ]1.0  + 

36.0  (T^/T^)Vs  exp  (-.15  - 0.5  T^/T^)]  cm/sec. 
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Here  consists  of  Spltzer  crossfield  resistivity  and  an  electron- 


neutral  resistivity  taken  from  Brown^  is  s value  obtained  from 

Liewer^^^;  is  a value  deduced  from  the  experiment  of  Zavoiskil, 

et  al.^®^  and  n is  from  the  parametric  model  of  Papadopoulos,^^^ 

'ee 

To  summarize  the  results,  by  choosing  the  initial  density  appro- 
priately, any  one  of  the  three  anomalous  resistivities  acting  alone  was 
capable  of  producing  the  desired  final  plasma  temperature  for  a beam 
current  density  in  the  range  3~5  kA/cm^.  (Classical  resistivity  alone 
required  an  order  of  magnitude  more  current  density.)  For  illustrative 
purposes,  we  choose  a case  in  which  the  electron-electron  counter- 
streaming  and  ion  acoustic  modes  are  present.  The  results  are  shown  in 
Fig,  5»  The  initial  electron  density  was  5^10  Icxsc  , the  neutral  density 
(H)  was  1.5xl0^^/cm^,  and  the  initial  temperatures  (T^(0)  = T^(0)  were 

3eV.  (The  results  are  insensitive  to  T (0),  but  do  depend  upon  n (0).) 

e e 

Figure  3a  shows  the  electron  and  ion  temperatures  as  a function  of  time 
during  and  shortly  after  the  beam  pulse.  The  initial  rapid  temperature 
rise  is  due  to  the  small  number  of  particles  sharing  the  large  energy 
input  (the  rf  s are  all  inversely  proportional  to  n)  at  early  time.  This 
temperature  rise  causes  rapid  ionization  (see  Fig.  3h)  which  slows  down 
the  heating  rate,  causing  the  temperature  to  drop.  Once  the  plasma  is 
> 90<^  ionized,  energy  input  again  goes  into  temperature  rise  Instead  of 
ionization,  but  now  at  a much  lower  rate.  The  heating  of  both  electrons 
and  ions  in  this  case  is  a result  of  assuming  the  energy  deposition  is 
partitioned  between  electrons  and  ions  as  given  by  Liewer*-  for  the 
two  instability  modes,  respectively.  Getting  full  ionization  depends  on 
the  beam  current  density,  as  shown  in  Fig.  k (the  initial  densities  and 
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Since  the  L/R  time  and  the  Implosion  time  (both  of  order  msecs)  are  long 
compared  with  the  magnetoacoustic  scale,  the  plasma  and  field  remain  In 
quasl-equlllbrlum  during  the  Implosion.  If  we  express  the  (azimuthal) 


current  density  in  terms  of  B by  Ampere's  law. 


f = - "oJ- 


and  rombin*^  this  with  (12)  and  (15),  we  find 


where  tj  represents  the  time-averaged  value  of  77*  Introducing  scaled 

variables  b = B/B  and  p = r/^wicb  A = have 

° \<TJ>  At/ 


g = - (<»>  / Tl)  ^ (l-b=)  . 


Here  (tj)  represents  the  spatial  average  of  tj. 

Depending  on  the  exact  microscopic  model  invoked  in  calculating 
the  anomalous  heating,  we  have 


T?  / <TJ>  = f [T,  B,  -g-  , etc.  ]. 


In  the  experiments  to  date,  the  induced  current  and  heated  plasma  are 
located  In  the  annular  region  originally  occupied  by  the  beam,  and  have 
radial  profiles  similar  to  the  density  profile  of  the  latter.  The 
observed'-  profiles  tend  to  be  broad,  occupying  a considerable  fraction 
of  the  total  system  radius.  It  seems  likely  that  the  plasma  profiles 
are  determined  mainly  by  that  of  the  beam,  not  by  the  form  of  tj.  One 


temperatures  were  held  equal  to  those  given  above,  and  the  same  anomalous 
resistivity  model  was  used  in  all  cases). 


j For  the  case  shown  In  Fig.  3,  the  plasma  current  dissipation, 

i followed  self-conslstently  In  time,  totals  10^  at  the  end  of  the  beam 

I 

pulse,  producing  500A/cm^  net  current.  This  Is  the  current  density  which 
persists  after  the  beam  leaves  the  plasma.  We  assume  che  1.^  MA  Gamble 
II  beam  Is  used  to  form  a 2^  cm  (average  radius)  annular  rotating  beam 

i 

with  a 10  cm  thickness  by  passing  through  a cusp.  If  V./V  ^ the 

0 z 

current  density  is  5 kA/cm^.  The  net  current  found  according  to  the 
calculation  described  above,  ^00  A/cm^,  Is  enough  to  reverse  the  3 kG 

(applied  field  on  axis  In  this  case.  This  configuration  has  a classical 

L/R  time  of  30  msec  for  T s 100  eV.  Since  the  current  density  Is  an 
^ order  of  magnitude  smaller  than  when  the  beam  was  present.  It  Is  reason- 

able to  expect  the  mlcrolnstabllltles  causing  anomalous  resistivity  will 
be  absent  during  this  time. 

The  beam-induced  return-current  heating  therefore  takes  place  prin- 
cipally over  the  lifetime  ^t  of  the  rotating  e-beam,  of  order  100  nsec. 
The  plasma  pressure  as  a function  of  position  Is  thus  given  by 


thus  concludes  that  rf  cannot  have  a very  strong  dependence  on  position. 
This  Is  In  contrast  with  the  situation  In  colllslonless-shock-heatlng 
theta  pinch  experiments  (with  and  without  bias  fields) In  which  the 
current  profile  Is  established  through  microscopic  processes  and  tends 
to  be  of  order  c/tOpj^  (Initially  1-2  orders  of  magnitude  smaller  than 
the  radius  for  the  parameters  we  are  considering). 

Thus  whenever  t]  is  large  enough  to  make  A comparable  with  the 
system  radius,  it  Is  reasonable  to  approximate  t)  = tj^  = constant,  whence 
from  (16)  b = - sin  p,  or 


B(r)  = - B sin 
o 


a) 


where  From  Eq,  (18)  it  follows  that  the  plasma  pressure  p(r) 

vanishes  at  r = r and  r = R^,  consistent  with  the  requirements  that  the 
density  vanish  on  axis  and  the  temperature  vanish  at  the  liner.  Figure  5 
shows  the  field  and  pressure  profiles  associated  with  Eq.  (18). 

Another  simple  solution  results  if  we  assume 


T,  pp  = 2^  (B^2  . g2) 


For  this  model,  Eq,  (16)  becomes 


4^  = • U j = const, 
dt 


Thus  the  current  density  is  uniform,  B varies  linearly. 


B = B (2r-R  )/R  , 

o o o 
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and  the  pressure  profile  Is  parabolic: 


p = f2B  2/u)  (r/R  ) (1-r/R  ).  (22) 

o ~ o o 

We  refer  to  these  two  models  as  the  sinusoidal  and  uniform  current  approx- 
imations. In  spite  of  the  difference  In  the  assumptions  on  the  functional 
form  of  T],  they  are  qualitatively  (and  quantitatively)  very  similar.  The 
numerical  results  described  In  the  next  section  depend  only  weakly  on 
which  Is  studied.  Again,  this  tends  to  justify  our  belief  that  the 
complicated  physics  embodied  In  the  definition  of  t)  Is  less  Important 
then  Its  average  magnitude.  This  Is  fortunate,  as  none  of  the  resisti- 
vities in  Eqs.  (8-11)  Is  well  approximated  by  a constant  or  by  Eq.  (I9). 

In  the  high  current  density  limit  for  which  Ionization  occurs 
rapidly  and  requires  only  a small  portion  of  the  total  energy  deposited 
by  resistive  heating,  the  electron  density  follows  the  Initial  neutral 
density  distribution.  With  uniform  Initial  gas  fill  the  Ion  acoustic  mode 
then  provides  a sinusoidal  current  distribution  self-conslstently. 
Electron-electron  two-stream  and  lower  hybrid  drift  Instabilities  provide 
current  distributions  that  are  Intermediate  between  the  sinusoidal 
and  uniform  current  models. 


III.  Adiabatic  Compression 

Let  the  system  be  compressed  by  reducing  the  liner  radius  from 

R to  R.  A differential  element  originally  at  r Is  displaced  to  a new 
o o 

position  r(r^),  where  0 ^ r ^ R,  and  undergoes  a volume  compression 
w(to): 
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■>»I.  >1  fn  I 


By  the  £rozen-ln  condition,  if  Pp  vanishes  at  r =0,  then 

o 

Pp(r^  = 0)  vanishes  at  all  subsequent  stages.  Furthermore,  r(0)  = 0 
by  symmetry.  We  can  parametrize  the  compressional  process  by  the 
volume  compression  \^o£  a plasma- free  differential  element  of  the 
system,  = w(0). 

The  instantaneous  radial  displacement  r of  a differential  element 
is  related  to  the  original  position  r by  Eq.  (25),  or  equivalently. 


dr^  _ 


= 2r 


j 2r 
rdr  - o 


o r dr 
o o 


Equations  (30)  and  (51),  together  with  the  initial  conditions  r(0)  = 0, 
w(0)  = w^,  can  be  integrated  numerically  to  provide  a parametric  definition 
of  the  function  w(r).  The  compressed-state  radius  follows  from 

R=r(R^).  (52) 

Instead  of  w^,  it  is  now  convenient  to  use  as  a measure  of  the  compression 
the  Inverse  scaled  system  radius, 

a = R^/R,  (55) 

defined  parametrically  by  (52). 

A knowledge  of  w(r)  = w [r(r^)l  suffices  to  generate  profiles  of 
density,  temperature,  pressure  and  magnetic  field  strength: 


n(r)  - n (r^)  w; 
o o 

T(r)  = T^(r^)  w^^^  ; 


I 

J 


V 


Figure  6 shows  how  6 and  the  average  beta  vary  with  a,  assuming  the 
uniform  current  model.  In  Fig.  7,  a log-log  plot  of  p vs.  a"^  is  shown. 
It  is  clear  that  the  compression  is  well  approximated  by  an  adiabatic 
law. 

2Y, 


^ eff 
pa  s const. 


(39) 


with  V ■ 1.82.  For  the  sinusoidal  model,  the  corresponding  value 


of  V is  1.85. 
eif 


In  the  course  of  the  implosion,  the  plasma  is  more  readily  com- 
pressed than  is  the  magnetic  flux.  As  a result,  the  "field-rich"  region 
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region  occupies  an  ever-lncreaslng  share  of  the  total  volume  at  the 
expense  of  the  "plasma-rich"  region.  This  is  reflected  in  the  marked 
reduction  in  € and  the  average  and  is  accompanied  by  an  Inward  shift 
in  the  relative  position  of  the  field-reversal  point. 

Assuming  a system  which  operates  with  a compression  value  a s 
we  find  a final  state  with  € = .50. 

For  this  state  the  radial  profiles  are  depicted  in  Figure  8. 

Assuming  uniform  Initial  density  and  writing  p = nkT,  we  can  calculate 
the  radial  dependence  of  n and  T separately,  as  shown.  Note  the  general 
steepening  of  all  gradients,  a feature  which  may  be  expected  to  reduce 
the  margin  of  MHD  stability  initially  present. 

Thus  far  we  have  worked  in  terms  of  scaled  variables.  If  we  specify 
the  final  conditions  according  to 

B = 1 MG  , 

max  * 


and  correspondingly. 


p = = 39  kbar, 

and  select  T to  be  in  the  range  of  thermonuclear  interest,  then  we 
max 

can  calculate  the  total  D-T  power  rate: 


P = 2.8x10"^'^  <^n^  t" 


-"4  V 

exp  (-19,914.1 

ave  ’ 


where  T is  temperature  in  keV  and  the  average  is  carried  out  over  the 
plasma  cross-section.  If  this  is  scaled  by  the  power  rate  which 

would  be  obtained  if  all  the  compresslonal  energy  resided  in  a plasma  of 
uniform  pressure  39  kbar  and  temperature  equal  to  and  the  results 
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plotted  against  temperature.  Fig.  9 Is  obtained.  The  curve  peaks  at 
approximately  I5  keV  (slightly  higher  than  the  maximum  of  the  unaveraged 
D>T  rate  for  constant  plasma  pressure),  with  values  of  order  10^  of  the 
optimum  ones. 

If  the  same  scaling  Is  used  with  the  average  of  the  bremsstrahlung 
power  rate, 

Pgj  = 5.3x10"^^  <n^T^  > ~ p^  <t"^>  , (41) 

Fig.  10  results.  Evidently,  for  a given  pressure  the  optimum  working 
temperature  lies  somewhere  above  1^  keV.  Lower  temperatures  (which  are 
naturally  easier  to  obtain)  penalize  the  system  In  terms  both  of  yield 
and  radiation  losses. 

IV.  Conclusions 

In  this  paper  we  have  discussed  (1)  the  form  of  the  reversed-fleld 
equilibrium  configuration  Induced  by  a rotating  e-beam  and  (2)  how  It 
changes  under  adiabatic  compression.  In  view  of  the  Idealizations 
employed,  our  conclusions  are  qualitative  and  somewhat  tentative. 
Nevertheless,  they  bear  closely  on  the  question  of  utilizing  the  reversed- 
fleld  theta-pinch  In  liner  reactor  designs,  and  so  are  of  practical 
Importance . 

First,  we  have  observed  In  Section  II  that  the  Initial  reversed- 
fleld  configuration  appears  to  be  determined  principally  by  the  pulse 
length,  current  density  and  radial  dimensions  of  the  relativistic  e-beam. 
(Confirmation  of  this  will  have  to  await  further  experimental  work.) 
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Our  heating  calculations  show  that  the  initial  configuration  does  not 
depend  strongly  on  the  details  of  the  heating  process,  i,e.,  the 
anomalous  resistivity  mechanism(s) . Further,  the  final  state  (for  a 
specified  g ) of  the  reversed-field  equilibrium  under  compression  is 
insensitive  to  the  exact  shape  of  the  initial  profiles  (Section  III) . 

The  second  conclusion  is  that,  without  optimization  of  the  initial 
configuration,  the  reversed-field  theta  pinch  leads  to  average  values 
of  beta  which,  though  substantial,  are  considerably  less  then  unity. 
Moreover,  because  of  the  nonlinear  functional  dependence  of  the  D-T 
power  rate,  the  yield  is  even  further  below  the  optimum  which  would  be 
achieved  if  only  uniform  plasma  were  compressed.  Against  this  must  be 
balanced  the  consideration  that  may  well  be  possible  to  create  sharper 
profiles  with  initial  £ higher  than  that  used  in  the  present  calcu- 
lations (.75)*  Further,  some  of  the  flux  at  large  radii  is  actually 
lost  by  diffusion  into  the  liner.  Although  such  losses  are  in  general 
undesirable,  it  is  important  not  to  count  the  energy  so  removed  twice- 
if  flux  is  lost  by  diffusion,  the  average  beta  must  increase  corre- 
spondingly. 

We  note  that  the  formalism  we  have  developed  is  appropriate  for 
calculating  the  adiabatic  radial  compression  of  an  arbitrary  reversed- 
field  theta  pinch.  It  may  also  be  used  directly  to  obtain  integrals 
(in  a linear  approximation)  of  total  D-T  yield,  radiation  losses, 
particle  and  heat  loss  and  magnetic  flux  diffusion,  as  in  the  example 
described  abv'>ve.  Evidently  there  is  no  need  to  restrict  consideration 
to  initial  states  for  which  an  analytic  representation  is  available. 

le 


The  present  formalism  will  be  applied  In  existing  codes  modeling  the 
self-consistent  dynamics  of  Imploding  liners.  We  have  shown  that  the 
pressure  of  such  a configuration  Is  well  described  by  an  ideal-gas 
adiabatic  law  with  an  effective  y 'ratio  of  specific  heats)  which  is 
weakly  dependent  on  the  shape.  The  successive  states  so  generated  can 
then  be  made  the  object  of  a study  of  both  hydromagnetic  and  micro- 
instabilities.  As  mentioned  previously,  both  types  depend  sensitively 
on  the  form  of  field  and  plasma  profiles.  If  the  initially  broad  field 
reversal  region  becomes  of  order  collective  phenomena  resembling 

those  in  collisionless  shocks  will  probably  occur,  giving  rise  to  an 
enhanced  resistivity  . 

In  order  that  the  profiles  calculated  according  to  the  assumption 
of  adiabatic  compression  not  be  smeared  out  by  diffusion,  we  must  have 

Dtj,  < A^,  (k-2) 

where  D is  the  radial  diffusion  rate,  t = a /V,  is  the  characteristic 

’ o o L 

d^^ell  time,  and  A is  the  sheath  thickness.  In  the  final  state  the  ions 

h 

are  strongly  magnetized  , where  is  the  ion  collision  time). 

Thus  we  can  write 

D = W , 

where  r^  = ^.5T^/B  cm  is  the  ion  gyroradius  (here  T is  in  keV, 

B in  kG) . Condition  (l)-2)  then  becomes 

V*  < -.C5  (A/a^)^  B^/T.  (1+5) 
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Taking  T = 10  keV,  B = 1 MG,  = 1 cm  and  = 3*10  cm/sec  yields 

i;*<  1.5  X 10®  (A/a^)^  sec"^.  (Ifl).) 

IfA/a^=0.1,  i;  < 1.5x10  sec  must  hold.  By  comparison,  the 

18 

classical  Ion  collision  rate  with  these  parameters  and  n = 2.5  x 10 
-5  5 -1 

cm  Is  1/  , « 9*10  sec  . Thus  If  the  collision  rates  in  the  final 

Cl 

plasma  state  are  classical,  our  model  should  be  adequate  to  provide  an 
indication  of  the  plasma/field  configuration. 
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Fig.  1 — Formation  of  a reversed-field-theta-pinch  by  a rotating  relativistic 
e-beam.  The  beam  enters  a drift  tube  filled  with  cool  gas  or  pleisma  (1),  re- 
versing the  field  on  axis.  The  beam-induced  return  current  decays  owing  to 
the  initial  high  plasma  resistivity,  heating  the  plasma  (2).  The  plasma  con- 
ductivity thus  increases,  so  that  when  the  beam  exits  (3),  a plasma  current 
is  induced  in  the  same  sense  as  that  of  the  beam.  This  persists  (4)  for  times 
long  compared  with  the  beam  lifetime,  since  the  plasma  has  become  hot  and 
highly  conducting. 


TIME  (ns) 


Fig.  3 — (a)  Electron  and  ion  temperature  variations  with  time  for  initial  electron 
density  n^,  = 3 X 10^®  cm~^,  neutral  H density  n„  = 1.5  X 10^®  cm~3  (20  m 
torr),  temperatures  Tj  = = keV,  current  density  J = 5 kA/cm^,  due  to  the 

action  of  classical,  electron-electron  two-stream  and  ion-acoustic  resistivities.  The 
nominal  beam  lifetime  is  90  nsec,  (b)  Electron  density  as  a function  of  time 
for  the  same  case. 
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u-u  RADIUS 

Fig.  6 — Initial  pressure  and  field  profiles,  (a)  Plasma  (solid  curve)  and  magnetic 
(broken  curve)  pressures  for  the  sinusoidal  case,  and  (b)  magnetic  field  strength 
for  the  same  case.  (Ck)ntinues) 


0-0  RRDIUS 

Fig.  5 (Ck>ntinued)  — Initial  pressure  and  field  profiles  (c)  plasma  and  magnetic 
pressures  for  the  uniform  current  case,  and  (d)  corresponding  field  profile. 


l.U 


COMPRESSION  RRTIO 


lb^.9 


Fig.  6 — Ratio  e of  plasma  energy  to  total  energy  (upper  trace)  and  average 
beta  (lower  trace)  as  a function  of  a,  the  radial  compression  factor,  for  the 
uniform  current  case 


u.u  RADIUS  U.uo 

Fig.  8 — Radial  profiles  for  the  uniform  current  case  with  a * 35  (e  =«  0.3), 
assuming  an  initially  uniform  density  profile  (cf.  Fig.  5):  (a)  plasma  (solid 
curve)  pressures  and  (b)  magnetic  field  (Continues) 


U.U  RADIUS 

Fig,  8 (Continued)  — Radial  profiles  for  the  uniform  current  case  with  a = 35 
(e  ^ 0.3),  assuming  an  initially  uniform  density  profile  (cf.  Fig.  5):  (c)  tempera- 
ture, and  (d)  density 
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Fig.  9 — Total  D-T  power  rate  for  the  state  shown  in  Fig.  8,  scaled  by  the  power 
which  would  be  obtained  if  all  the  compressional  energy  had  gone  its  plasma 
heating.  The  values  shown  are  for  temperature  T = 5,  10,  ...,  25  kev,  assuming 
constant  pressure  and  Maxwellian  ion  distributions. 


